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IDENTIFICATION AND USE OF INFORMATIVE SEQUENCES 

CROSS REFERENCE TO RELATED APPLICATIONS 

[0001] The present application claims the priority benefit of pending U.S. Provisional Patent Application 
No. 60/441,745 filed January 23, 2003 and pending U.S. Provisional Patent Application No. 60/441,806 
filed January 23, 2003. The complete disclosure of both applications is hereby incorporated herein by 
reference. The present application also hereby incorporates herein by reference the entire disclosure of 

pending utility patent application 10/ , titled Method and System for Identifying Biological Entities 

in Biological and Environmental Samples, filed January 23, 2004. 

STATEMENT REGARDING FEDERALLY SPONSORED RESEARCH OR DEVELOPMENT 

[0002] Embodiments of the invention may have been conceived or first actually reduced to practice in 

the performance of work under the following Government contracts: . As a 

result, the Government may have certain rights in those inventions. 

BACKGROUND 

[0003] Field of the Invention . Embodiments of the invention relate to the identification of genomic 
sequences that are informative of the biological characteristics (e.g., presence, abundance, virulence, 
genetic modification) of a sample, along with systems and methods of using such sequences for gathering 
information on one or more sets of organisms present in the sample. Specific embodiments relate to 
microbial organisms. 

[0004] Description of Related Art . Genes, natural units of hereditary material, are the physical basis for 
the transmission of the characteristics of biological entities from one generation to another. The basic 
genetic material is fundamentally the same in all biological entities. It consists of chain-like molecules of 
nucleic acids (deoxyribonucleic acid (DNA) in most organisms and ribonucleic acid (RNA) in certain 
viruses) and is usually associated in a linear or circular arrangement that, in part, constitutes 
chromosomes and extra-chromosomal elements, such as micro-chromosomal bodies. The entire 
hereditary material in a cell is called the "genome." In addition to the DNA contained in the nucleus, an 
organism's cells contain DNA in other locations within those cells, e.g., bacteria also contain some DNA 
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in plasmids, plants also contain some DNA in plastids, animals also contain some DNA in mitochondria. 
A set of biological entities, such as a species, has a genome, e.g., the complete sequence of genes 
characteristic of the set. Some portions of the genome are unique to the particular set, e.g., set-unique 
sequences. Example sets include strain, species, genus, family, group, clade, and other ad hoc sets. 

[0005] Historically, the theory, principles, and process of classifying biological entities into sets {e.g.,, 
taxonomic classification) is based on the work of seventeenth century biologist Carl Linnaeus. Linnaeus 
created the taxonomy system of kingdom, phylum, class, order, family, genus, and species. Known as the 
Linnean system, this rank-based taxonomy is still in use today. Other basis for classifying organisms 
have been proposed, including some based on phylogeny, i.e., the evolutionary development of biological 
entities. For example, in contrast to the rank-based codes , the PhyloCode will provide rules for the 
express purpose of naming clades and species through explicit reference to phylogeny. See e.g., 
http://www.ohiou.edu/phylocode/index.html, accessed January 14, 2004. 

[0006] As noted by the National Center for Biotechnology Information (NCBI) at 
htpp://www.ncbi.nlm.nih.gov/BLAST/blast_overview.shtml (accessed January 5, 2004), BLAST ® 
(Basic Local Alignment Search Tool) is a set of similarity search programs designed to explore genetic 
sequence databases available through NCBI. The BLAST programs have been designed for speed, with a 
minimal sacrifice of sensitivity to distant sequence relationships. The scores assigned in a BLAST search 
have a well-defined statistical interpretation, making real matches easier to distinguish from random 
background hits. BLAST uses a heuristic algorithm that seeks local as opposed to global alignments and 
is therefore able to detect relationships among sequences that share only isolated regions of similarity. 
The Expected Value ("is") as noted in BLAST search results is a parameter that describes the number of 
hits of the type shown that one can expect to see just by chance when searching a database of a particular 
size. It decreases exponentially with the Score ("5") that is assigned to a match between two sequences. 
E can be interpreted as the random background noise that exists for matches between sequences. For 
example, an E value of 1 assigned to a hit can be interpreted as meaning that in a database of the current 
size, one might expect to see one match with a similar score simply by chance. This can be interpreted to 
means that the lower the is-value, or the closer it is to "0", the more significant the match is. 
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[0007] Embodiments of the invention include systems and methods to identify of genomic sequences 
that are informative of the biological characteristics (e.g., presence, abundance, virulence, genetic 
modification) of a sample, along with systems and methods of using such sequences for gathering 
information on one or more sets of organisms present in the sample. Specific embodiments relate to 
microbial organisms. 

[0008] A method for identifying genomic sequences unique to a set of organisms includes obtaining 
genomic data characteristic of the set. The sequences are formatted into query-length sequences; each 
query-length sequence being of a format compatible with a similarity search engine such as BLAST. A 
selected genomic database, such as those maintained by NCBI (or a more restricted database)is searched 
using the query and the similarity search engine. The results of the search are parsed for those sequences 
showing uniqueness to the set. 

[0009] The invention includes computer programs for identifying genomic sequences unique to a set of 
organisms. These programs can be carried one or more computer-readable media and include a genomic 
data interface module. The genomic data interface module is operable to couple to a source of genomic 
data to receive genomic data characteristic of the set. A formatting module formats the received genomic 
data into query-length sequences, where each query-length sequence is formatted compatible with a 
similarity search engine such as BLAST. A search interface module interfaces with the similarity search 
engine to submit the query-length sequence to the search engine. A search results parsing module parses 
the results of the search for those sequences showing uniqueness to the set. 

[0010] Embodiments of the invention include methods for identifying oligonucleotide sequences unique 
to a set of organisms. In those methods, after identifying unique genomic sequences (typically, but not 
necessarily, as described above), unique genomic sequences are divided into target-length oligonucleotide 
sequences; preferable a number of sequences starting with the beginning of the genomic sequence and 
progressing one nucleotide at a time. The oligonucleotide sequences are properly formatted and searched 
to a selected database using a similarity search engine such as BLAST. The results are parsed for those 
oligonucleotides showing uniqueness to the set. A computer program product can implement this method 
in modules similar to those above. 



-3- 



PATENT SAIC0086-US 

[0011] The invention includes methods for inferring both unique genomic sequences and unique 
oligonucleotide sequences. In those methods, genomic data characteristic of a first set of organisms is 
obtained. The data is formatted into at least one query-length sequence, each query-length sequence 
being of a format compatible with a similarity search engine. A selected genomic database, e.g., an NCBI 
database, is searhed using the query and the similarity search engine such as BLAST. The results are 
parsed for those sequences not associated with the first set of organisms, but showing similarity beyond a 
threshold. 

BRIEF DESCRIPTION OF THE DRAWINGS 

[0012] Figure 1 is a flowchart describing, in conjunction with portions of the written description, 
methods of the present invention. 

[0013] Figure 2 is a notional illustration of a target oligonucleotide window in the context of a unique 
genomic sequence of the present invention. 

[0014] Figure 3 is a notional illustration of differential hybridization of several unique oligonucleotides 
of the present invention between two strains of E. coli. 

[0015] Figure 4 is an illustration of a decision tree of the present invention. 

DETAILED DESCRIPTION OF THE INVENTION 

[0016] As required, detailed embodiments of the present invention are disclosed herein. However, it is 
to be understood that the disclosed embodiments are merely exemplary of the invention that may be 
embodied in various and alternative forms. The figures are not necessarily to scale, and some features 
may be exaggerated or minimized to show details of particular components. Therefore, specific structural 
and functional details disclosed herein are not to be interpreted as limiting, but merely as a basis for the 
claims and as a representative basis for teaching one skilled in the art to variously employ the present 
invention. 

[0017] For purposes of the invention disclosure, the term "primer" includes a short pre-existing 
polynucleotide chain to which new nucleotides can be added by DNA or RNA polymerase. 
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[0018] The term "randomly amplifying" includes increasing the copy number of a segment of nucleic 
acid in vitro using random primers. 



[0019] "Amplicon" refers to DNA that has been manufactured utilizing a polymerase chain reaction 
(PCR) where a set of single stranded primers is used to direct the amplification of a single species of 
DNA. 

[0020] "Biological entity" describes a biological element, cellular component , or organism that exists as 
a particular and discrete unit. This includes, but is not limited to gene, transgene, oncogene, allele, 
protein, DNA, RNA, mitochondria, pathogenic trait, vector, plasmid, clone, Acytota, prokaryotes, 
eukaryotes, Protista, Fungi, Plantae, Animalia and Monera, or any mixture thereof. For simplicity, 
"organism" is used interchangeably herein with "biological entity," and includes sub-organism entities. 
"Set of organisms" includes sets of one biological entity. 

[0021] A "sample" may be from any source, and can be a gas, a fluid, a solid, a biological sample, an 
environmental sample, or any mixture thereof. 

[0022] "Nucleic acids" means RNA and/or DNA, and may include unnatural bases. 

[0023] The "unique oligonucleotide sequence" generally identifies a nucleic acid sequence for which the 
sequence is known and determined unique to a set of organisms. In some embodiments, unique 
oligonucleotide sequences are more than 30 nucleotides in length. 

[0024] The terms "unique genomic sequence" and "unique sequence" are interchangeable in the 
invention and refer generally to a sequence of nucleic acids that are specific to a set of organisms. 

[0025] In the literature there exist at least two confusing nomenclature systems for referring to 
hybridization partners. Both use common terms: "probes" and "targets." For the purpose of this 
disclosure, a "target" is the known oligonucleotide sequence (preferably set-unique), whereas a "probe" is 
the nucleic acid sample whose characteristic(s) (e.g., identity, abundance, virulence) is being detected. 
"Probe" includes any single stranded nucleic acid sequence, molecule, genomic sequence, or amplicon 
which typically is labeled. Probes can hybridize to a target if sufficient complementarities exist. Note 
that labeling can be implemented at various stages in either the probe or target or both, as known to those 
skilled in the art. 
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[0026] The terms "microarray" and "array" are interchangeable and include a set of miniaturized 
chemical or biological reaction areas that may also be used to test DNA, DNA fragments, RNA, 
antibodies, or proteins. 

[0027] A "labeled" or "detectable" nucleic acid is a nucleic acid that can be detected. The term 
"detection" refers to a method where analysis or viewing of the detectable nucleic acid is possible visually 
or with the aid of a device, including, but not limited to microscopes, fluorescent activated cell sorter 
(FACS) devices, spectrophotometers, scintillation counters, and fluorometers, devices using mass 
spectrometry, devices using radio isotopes. 

[0028] "Hybridized" means having formed a sufficient number of base pairs to form a nucleic acid that is 
at least partly double-stranded under the conditions of detection. The term "hybridization" refers to the 
process by which two complementary strands of nucleic acids combine to form double-stranded 
molecules. 

[0029] The term "complementarity" refers to a property conferred by the base sequence of a single 
strand of DNA or RNA that may form a hybrid or double stranded DNA:DNA, RNA:RNA or DNA:RNA 
through hydrogen bonding between base pairs on the respective strands. Adenine (A) usually 
complements thymine (T) or uracil (U), while guanine (G) usually complements cytosine (C). 

[0030] The identification and characterization of microorganisms present in the environment has 
historically been accomplished by exploiting a variety of biological, immunological, biochemical, and 
genetic differences between organisms. Nucleic acid-based diagnostic methods have been developed that 
are specific for a single organism or small sets of organisms. PCR-based assays are typically performed 
by designing oligonucleotide primers that amplify organism-specific fragments of DNA. These 
fragments are subsequently detected by methods such as gel-electrophoresis, real-time PCR, or 
hybridization to either a membrane or microarray. A limitation of these existing assays is that although a 
positive result is informative for a specific organism or organism set, a negative result typically provides 
little or no information about the organism(s) under investigation. 

[0031] Though it is possible to multiplex primers for numerous amplifications in a PCR for the 
concurrent identification of a variety of organisms, it is non-trivial to design compatible multiple primer 
pair sets that function in a single amplification reaction. Thus, the number of microorganisms that can be 
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detected or otherwise characterized concurrently with this type of multiplex reaction is relatively small. 
Techniques such as real-time PCR and quantitative PCR are limited in the number of primer sets that can 
be used in a single amplification reaction and in the number of fluorescent molecules available for 
labeling DNA molecules and detection. 

[0032] In a method for detecting and distinguishing between various species and strains of viruses, viral 
RNA is reverse transcribed from semi-random primers, amplified by specific primers and then labeled 
with fluorescent nucleotides in a non-amplifying reaction. The labeled nucleic acids are then hybridized 
to microarrays that have been spotted with virus and strain-specific oligonucleotides that represent the 
entire genomes of these organisms. The resulting hybridization pattern discriminates between viruses 
represented on the array. However this approach is not directly translatable to fungi and bacteria. The 
relatively large size (3-5 million bases) and complexity of bacterial and fungal genomes, as compared to 
most viral genomes, represents an obstacle in the ability to identify oligonucleotides that are species and 
strain specific. In addition, it is not feasible to synthesize and spot every possible oligonucleotide for 
every microbial genome onto a microarray. 

[0033] Bioinformatic tools such as BLAST, are intended to identify similarities between sequences. 
While similarities between the sequences of organisms are useful in some types of analysis, the 
differences between genomes can be useful in the identification and characterization of microorganisms. 
Unfortunately, bacterial and fungal genomes are so large that it is resource-intensive to subtract common 
sequences in order to identify unique sequences from all known genomes. Frequently only small 
fragments of genomic sequences have been identified as unique are available for identification of an 
organism. The increasing number of genomes that have been, or will soon be, sequenced is one incentive 
for identifying large fragments of known genomes. 

[0034] Current nucleotide-based methods of identifying organisms typically rely on primer-requiring 
multiples PCR methods, or oligonucleotide microarrays that utilize the limited amount of ribosomal genes 
(approximately 1% of the genome), or costly shotgun approaches directed at entire genomes. Ribosomal 
DNA is highly conserved across higher level sets - therefore the amount of unique sequence is limited. 

[0035] Approximately 300 microbial genomes have been partially or completely sequenced to date. In 
spite of this wealth of information, existing methods of detection and characterization of microbes are 
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limited by the availability of unique sequence information from the genomes of these organisms. 
Typically, only small fragments of genomic sequences are identified as unique. 

[0036] Current DNA amplification approaches to identify microorganisms are limited in terms of the 
number of sequences that can be identified concurrently. In vitro, two methods are used to multiplex, or 
identify multiple sequences concurrently. Both are limited by the challenge of generating specific primer 
pair sets that work well together in a single reaction mixture. One method for assessing which amplicons 
are produced in a multiplex PCR is to run the amplification product on a gel and to discriminate the 
various amplicons based on molecular size. The number of size indicative bands that can be resolved on 
the gel is a limiting factor for this approach. 

[0037] Another method, real time PCR, uses different fluorescent tags to identify specific amplicons in a 
multiplex PCR. The number of amplicons that can be resolved using this approach is limited by the 
number of different fluorescent tags available for probes used in the reaction. Current limitations for 
fluorescent resonance energy transfer methods, such a Taqman and molecular beacons are about four 
amplicons for a single reaction. Thus the limitations are at least two-fold: the first is a compatibility issue 
regarding the use of multiple sets of unique primers; and the second is resolution of the amplified 
products. 

[0038] Unique genomic sequences in a set of organisms (which set can consist of one or more biological 
entities) may include both coding and non-coding sequences. Coding sequences are sequences that are 
further processed into proteins or polypeptides, typically performing a single function. These sequences 
are frequently conserved across genus and species. Conserved coding sequences can include genes that 
code for enzymatic elements, structural elements, virulence factors or developmental specific functions 
and processes. 

[0039] Non-coding sequences are sequences that are not further processed and do not appear to possess a 
known function at this time. These sequences may be contained in a portion of the genome that contains 
unique coding sequences as well as between conserved coding sequences. Since non-coding sequences 
do not provide a known function, they are frequently overlooked as unimportant genomic material. These 
unique non-coding sequences can be used to identify an organism, just as unique coding sequences are 
used. Informative sequences can reflect a variety of features e.g. structural, functional, metabolic, 
virulence. Set-unique sequences can be coding or non-coding sequences. Set-unique sequences (coding 
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or non-coding) can be inferred (see below) or found by searching through fully sequenced genomes. 
Partially sequenced genomes typically focus on coding sequences. Combinations of sequences that are 
not necessarily individually set-unique can also be informative. Sequences unique to sets above the 
species or strain level can be bio-informative, e.g., used in analyzing a sample for information about the 
organisms in the sample as described below. 

[0040] Embodiments of the present invention include methods and systems for the identification of 
genomic sequences that are informative of the biological characteristics (e.g., presence, abundance, 
virulence, genetic modification) of a sample. Referring to Figure 1, a illustrative method 100 of the 
present invention is shown. 

[0041] In the illustrated embodiment, a subset of the genomic data 105 of the organism under 
investigation is obtained. The subset 105 can be obtained from known genomic data source 10, e.g., 
UniGene, GenBank, European Molecular Biological Laboratory (EMBL), among other sources. 
Genomic data can also be obtained as sequence derived from in vivo or in vitro experiments 20 such as 
PCR and enzymatic digestion. A preferred subset of genomic data is the entire genomic sequence of an 
organism. 

[0042] In some embodiments, the obtained genomic data is preprocessed 110. Each aspect of 
preprocessing can be performed as needed or desired. 

[0043] In preferred embodiments, if necessary, the genomic data subset is converted from its native 
format, e.g., standard GenBank annotated format, to a format compatible with subsequent steps. In some 
embodiments, where GenBank annotated form at is used, the genomic data is converted to FASTA format 
to support a subsequent BLAST search. 

[0044] Preprocessing can involve removing or masking portions of the genomic data that are judged not 
likely to have informative value. In preferred embodiments, these portions are removed, partly to make 
the subsequent search more efficient. This can include sequences known to be conserved with respect to 
the organism set under investigation (though some conserved sequences can be bioinformative, and 
sequences conserved at one tier in a taxonomy may be unique in another tier), repeats, inverted repeats, 
long terminal repeats, sequences otherwise known to be not favorable for hybridization. 
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[0045] In preferred embodiments, genomic data is divided into query-length sequences 115. Some 
embodiments start with sequences of 1000 bases in length. The speed of the search is one factor in 
selecting the size of the initial query sequence. Subsequent iterations, described below, divide the initial 
query length further. 

[0046] In further preprocessing, query-length genomic sequences were realigned with the genome from 
which they were generated in order to determine the start and stop point of each query length sequence 
within the genome. Any annotations within the genome in the region containing the query length 
genomic sequence were transferred to the query length genomic sequence. Annotated regions include 
sequences known to have a specific biological function such as protein coding regions, biologically active 
RNA encoding regions, promoter and regulatory elements, spacing elements within operons, protein 
binding sites, and the like. 

[0047] Note that if an entire genome was obtained, and no preprocessing performed, the query-length 
sequence 115 is the entire genome of the organism under investigation. In some embodiments, all the 
genomic data available for the organism under investigation is obtained, all preprocessing steps are 
completed, resulting in annotated query-length sequences of 1000 bases that do not include conserved 
sequences, repeats of various types, or sequences having characteristics that otherwise make them not 
amenable to subsequent steps. 

[0048] In preferred embodiments, the query length sequence (preprocessed or not) is used as a query to a 
similarity search program 120, e.g., BLAST. The query is directed to a selected database 125 of genome 
data. In some embodiments, the selected database is limited to organisms of the same general type as that 
under question, in order to increase search efficiency over what it would be were the search directed to a 
full database containing a broader variety of organisms. In some embodiments, the query is directed to 
the NCBI nr database. For example, if only microbial organisms were under investigation, the selected 
database 125 would be a database of microbial genomic data - broader databases including, for example, 
mammalian genomic data, would be avoided at this stage. In these circumstances, a subsequent search 
against the broader database is preferred in order to confirm the set-uniqueness of these initial results. In 
some embodiments, query-length sequence is removed from the selected database, while in other 
embodiments, results showing homology to the query itself are either ignored, or taken as confirmation of 
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the validity of the query with respect to the organism under investigation. Table 1 lists an exemplary 
query for a strain of Escherichia coli. 



>NC_000913_2 9__part354 Escherichia coli K12, complete genome. 

TATGTTTTAAATAACTAATTGGTCGGGTTAGTGCATCCGGCTTTCTTTATATTCGCCAGA 

AGGATTTATTATGCAAAGGAAAACTCTATTGTCGGCCTGTATTGCATTAGCTCTGAGTGG 

TCAGGGTTGGGCGGCAGATATCACAGAGGTAGAAACCACCACAGGTGAAAAGAAAAATAC 

CAATGTGACTTGTCCGGCAGACCCAGGAAAACTCAGTCCGGAAGAGCTTAAACGCTTACC 

CTCTGAATGCTCTCCTTTAGTCGAACAAAACCTGATGCCATGGCTTTCCACAGGCGCTGC 

TGCGTTAATCACGGCCTTAGCCGTAGTGGAACTAAACGACGATGATGATCATCATCATCG 

CAACAATTCTCCACTCCCACCGACACCCCCTGATGATGAATCAGACGACACTCCAGTTCC 

CCCAACTCCTGGCGGAGATGAGATAATACCGGACGATCCGGATGATACGCCTACACCTCC 

CAAACCGGTTTCGTTTAATAATGACGTTATTCTCGATAAAACAGAAAAAACGTTAACTATT 

CGCGATTCAGTTTTTACTTATACCGAGAATGCTGACGGGACTATATCTCTGCAAGATAGCA 

ATGGTCGTAAGGCAACGATTAATCTTTGGCAGATTGATGAAGCGAATAACACTGTTGCCCT 

TGAAGGGGTGAGCGCAGATGGCGCAACGAAGTGGCAATATAATCACAACGGTGAGCTTGTT 

ATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATT 

CCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGT 

CAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGC 

ACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGT 

CAATAACGAAGGCGAGAGCACCATCACCAAC 



TABLE 1 



[0049] Preferred embodiments parse 130 the similarity search program output 125 to identify sequences 
lacking significant similarity with other organisms in the selected database, e.g., unique genomic 
sequences 132. This is counter to the typical use of such search programs. For the purpose of this 
disclosure, "unique" or "uniqueness" is a function of thresholds, preferably controlled by the user, 
regarding identity, homology, score, expected (E) value and the length of the unique sequence under 
consideration. Identity, score, expected value are data returned in a typical BLAST search. In some 
embodiments, lacking significant similarity, e.g., "unique," means no BLAST hits or hits with a E-value 
less than le-5. 

[0050] In practice, computational resources are finite, so the selected database may range from a 
database of all fully or partially known genomes to a narrower database such as known microbial 
genomes. Directing the initial query to a database of less than all available genome data, while 
computationally economical, can make it advisable to BLAST the candidate sequences (e.g., in preferred 
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embodiments, those genetic sequence segments found to be unique) against the broader databases, e.g., 
the NCBI nr database, to detect homology with known genomes. 

[0051] Table 2 is an extract from an exemplary BLAST output for the sequence in Table 1. Note that 
significant alignments with E values less that le-47 are all E. coli. This confirms that the query sequence 
is sufficiently "unique" with respect to E. coli to be informative. 



BLASTN 2.2.7 [Jan- 02 -2 004 ] 

RID: 1074618061-23468-21421924927. BLASTQ4 
Query= (1000 letters) 



Database: All GenBank+EMBL+DDBJ+PDB sequences (but no EST, STS, 

GSS, or phase 0, 1 or 2 HTGS sequences) 2,017,250 sequences; 9,771,119,756 

total letters 

Distribution of 64 Blast Hits on the Query Sequence 



Sequences producing significant alignments: 

gi | 1787665 |gb|AE000237 , 1 |AE000237 Escherichia coli K12 MG16 . 
gi j 1742273 | dbj |D90778 . 1 | E.coli genomic DNA, Kohara clone #. 
gi|2610794l|gb|AE016760.l| Escherichia coli CFT073 section . 
gi j 41829 | emb |X62680 . 1 | ECIS2IS30 E.coli insertion sequences . 
gi|24431548|gb|AC125613.2| Homo sapiens 3 BAC RP11-379M20 (. 
gi | 18159853 |gb|AE009803 . 1 |AE009803 Pyrobaculum aerophilum s. 
gi j 6067150 | gb | AC008015 . 5 | AC008015 Homo sapiens chromosome 1. 
gi | 21644686 | dbj | AP004363 . 3 | Oryza sativa (japonica cultivar. 
gi | 10862756 | embj AL354852 . 12 | Human DNA sequence from clone . 



Score 


E 


bits) 


Value 


1982 


0.0 


1982 


0.0 


1191 


0.0 


198 


le-47 


44 


0.50 


42 


2.0 


42 


2 . 0 


42 


2.0 


42 


2.0 



TABLE 2 



[0052] At this point, those sequences (less than or equal to query-length) that show homology/identity 
with other organisms below a threshold can be identified as unique to the set for which they were 
searched. Such sequences have utility to confirm the presence of at least one member of the set, 
primarily, but not exclusively in a Bioinformatic setting. Sequences unique to higher level sets are 
identified by searching for commonalities between sequences within a classification. These common 
sequences are then searched against the appropriate database. Such sequences, with optional annotation, 
can be realigned with the genomic data, annotated with the information gained saved to a database of 
unique genomic sequences 132, or added to the growing knowledge base of the genome of the organism 
under investigation. 



- 12- 



PATENT SAIC0086-US 

[0053] In some embodiments, where a certain minimum number of total unique bases , e.g., 10,000, is 
desired across all unique genomic sequences, and less than this number are determined to be unique, the 
genomic data can be preprocessed again, this time dividing the genomic data into smaller query-length 
sequences. The smaller query length sequences are then searched against the target database. Control 
line 133 indicates this path. 

[0054] The output of the similarity search program can also be used to identify further query-length 
sequences or candidate sequences for organism(s) other than the organism(s) under investigation. For 
example a first query-length sequence may show high homology/identity only against the particular strain 
it was derived from. But the query sequence might also show some homology to a related strain. Such 
sequences can be referred to as inferred sequences 134. The portion of the related strain where limited 
homology is detected can be searched 120 as a query-length genomic sequence 115 (by being searched 
against the selected database 125) to confirm its identity as a unique genomic sequence 132 for the related 
organism(s); or it can be treated as a candidate sequence for evaluation (discussed below) where target- 
length oligonucleotides are evaluated for amenability to hybridization. Exemplary inferred sequences 
have sufficient homology to the related sequence to be indicated by a BLAST search, but not sufficient to 
cross-hybridize with oligonucleotides derived from the related sequence. 

[0055] Referring to Table 3, a search against the NCBI nr database, using as a query a Vaccinia 
sequence found to be unique by a method of the present invention, identified two candidate sequences 
A01, A02 (a Vaccina strain and the complete Vaccinia genome) with 100% identity over the entire query 
sequence; five non-Vaccinia sequences A03 - A07 with identity ranging from 96% to 92% over portions 
of the query sequence; one non-Vaccinia sequence A08 with 100 identity over a small portions of the 
query sequence; and at least seven non-microbial sequences A09 - A15 with E values greater than three 
for short portions of the query sequence. 

[0056] The first group confirms that the query sequence is part of both the Vaccinia strain and complete 
genome. The second and third groups identify sets of organisms with significant homology, e.g., E value 
less than le-5, to the set-unique Vaccinia sequence. Preferred embodiments of the invention infer that the 
second and third group of sequences come from unique regions of the genome of those organism sets. 
Such inferred sequences preferably undergo evaluation and validation as described herein. 
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BLASTN 2.2.4; RID: 1036169670-05727-22152 






Query= (160 letters) 






Database: All GenBank+EMBL+DDBJ+PDB sequences (but no EST, STS, GSS, or phase 0, 1 or 2 
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gipjjji / 1 gD|ivi j jKjl /. 1 1 v /w^v^o vaccinia virus, complete genome 


117 


zc-oh 


gipuyoyoz|emu| i i io*+z.i|L^ vvxtviyu v^owpoA virus sixain oivi-!/v/ 




Sp-70 


gipuy /ui j|emu| Y 1 jUjj. i|uv y i jU3D v^owpox virus strain LrKi-yu... 


77H 
Z /U 


De-/u 


gi|20152989|gb|AF482758.1| Cowpox virus strain Brighton Red... 


252 


le-64 


gij 1 84829 1 3 |gb| AF43 8 1 65 . 1 1 Camelpox virus M-96 from Kazakhs... 


228 


2e-57 


gijl 97 17929|gb| AY009089.1| Camelpox virus CMS, complete genome 


220 


4e-55 


gi|22123748|gb AF012825.2| Ectromelia virus strain Moscow, ... 


80 


9e-13 


gi|14574206|gb|U23449.2| Caenorhabditis elegans cosmid K06A... 


38 


3.2 


gi|687828|gb|U21318.1| Caenorhabditis elegans cosmid K03H9,... 


38 


3.2 


gi|12000447|gb|AC084754.14| Homo sapiens 12p BAC RP11-874G1... 


38 


3.2 


gi|17534934|ref|NMJ)62895.1| Cuticulin precursor 


38 


3.2 


gi|18250549|emb|AL627429.8| Human DNA sequence from clone R... 


38 


3.2 


gij 16973060|emb|AL590 101.9| Human DNA sequence from clone R... 


38 


3.2 


gi|23337297|emb|AL7323 17.13| Mouse DNA sequence from clone ... 


38 


3.2 


TABLE 3 1 



[0057] Unique and inferred unique genomic sequences can be identified using the method described 
herein for a number of other biological entities including, but not limited to; Anthrax (Bacillus anthracis), 
Botulism (Clostridium botulinum toxin), Brucellosis (Brucella species), Burkholderia mallei (glanders), 
Burkholderia pseudomallei (melioidosis), Chlamydia psittaci (psittacosis), Cholera (Vibrio cholerae), 
Clostridium perfringens (Epsilon toxin), Coxiella burnetii (Q fever), E. coli 0157:H7 (Escherichia coli), 
Emerging infectious diseases such as Nipah virus and hantavirus, Food safety threats (e.g., Salmonella 
species, Escherichia coli 0157:H7, Shigella), Francisella tularensis (tularemia), Ricin toxin from Ricinus 
communis (castor beans), Rickettsia prowazekii (typhus fever), Salmonella Typhi (typhoid fever), 
Salmonellosis (Salmonella species), Smallpox (variola major), Staphylococcal enterotoxin B, Variola 
major (smallpox), Viral encephalitis (alphaviruses e.g., Venezuelan equine encephalitis, eastern equine 
encephalitis, western equine encephalitis), Viral hemorrhagic fevers (filoviruses e.g., Ebola, Marburg and 
arenaviruses e.g., Lassa, Machupo), Yersinia pestis (plague). 
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[0058] Referring again to Figure 1, unique genomic sequences were realigned 130 with the genome from 
which they were generated in order to determine the start and stop point of each fragment within the 
genome. In preferred embodiments, annotations within the genome in the region containing the unique 
sequence was transferred to the unique sequence. Annotated regions include sequences known to have a 
specific biological function such as protein coding regions, biologically active RNA encoding regions, 
promoter and regulatory elements, spacing elements within operons, protein binding sites, etc. 

[0059] In some embodiments of the present invention, the process of obtaining genomic data, 
preprocessing the data, querying the selected database(s) and parsing results to identify candidate 
genomic sequences is implemented as a computer program product. In these embodiments, a plurality of 
organisms and sets of organisms can be investigated concurrently. Computer program products of this 
invention include the ability to indicate the organism(s)/set of organisms of interest, indicate the selected 
database, set thresholds for identifying inferred sequences, direct the handling for inferred sequences, set 
thresholds for identifying unique sequences, direct the handling for unique sequences, output unique 
sequences for evaluation for oligonucleotides amenable to hybridization. Intermediate and final results 
can be made available for user inspection. In preferred embodiments, such computer program products 
are in network communication, e.g., via the Internet with selected databases and databases of genomic 
data such as those available through NCBI via http://www.ncbi.nlm.nih.gov. Operator interface to such 
computer program products is preferably provided through graphical user interface (GUI) technologies as 
known to those skilled in the art. Such computer program products can be configured to operate in a 
centralized fashion or may be distributed over platforms on a network. Some embodiments are 
configured as an application service provider (ASP) accessible by network devices through a browser 

[0060] Both unique genomic sequences 132 and select inferred unique sequences 136 are evaluated 140 
for subsets e.g., target-length oligonucleotides, that are amenable to hybridization. The evaluation is done 
in a target-length oligonucleotide window derived from the query length sequence, and preferably moved 
one base at a time through the query-length sequence. Figure 2 presents a notional representation 200 of 
the incremental progression of a target-length oligonucleotide window 235 through a unique genomic 
sequence 232 one nucletotide 236 at a time. In preferred embodiments, the increment of progression is 
one nucleotide at a time, but in other embodiments different constant or variable progressions can be used 
to take advantage of the inherent properties of the unique genomic sequence. Target-length 
oligonucleotides are evaluated for, among other characteristics, GC content, metlting temperature (TJ, 
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repetitive elements, availability of primer amplification sites, avoiding secondary structures such as 
hairpins and duplexes. In some embodiments this functionality is provided using a program such as 
OLIGO 6 (Molecular Biology Insights, Inc., Cascade CO). In other embodiments, this functionality is 
incorporated into a computer program product of the invention. 

[0061] OLIGO is a multi-functional program that searches for and selects oligonucleotides from a 
sequence file for polymerase chain reaction (PCR), DNA sequencing, site-directed mutagenesis, and 
various hybridization applications. It calculates hybridization temperature and secondary structure of 
oligonucleotides based on the nearest neighbor thermodynamic values. It is also a good tool for 
construction of synthetic genes, finding an appropriate sequencing primer among those already 
synthesized, finding and multiplexing consensus primers and probes, and even finding potential 
restriction sites in a protein. 

[0062] In preferred embodiments, oligonucleotides of approximately 50 bases in length are derived from 
the candidate sequences. A preferred range for oligonucletide lengths is 25- 100 with a range of 50 - 70 
being more preferred. Factors that go into determining a range and a preferred value include the ability to 
synthesize the oligonucletoide, the desired hybridization temperature of the microarray, balancing melting 
temperature of the various oligonucleotide against the GC content of the molecule and the possible 
chemical composition of the hybridization solution used on the microarray. As these factors change, the 
preferred length of oligonucleotides will also change. In some embodiments, target-length 
oligonucleotides are chosen based on their melting temperature T m of 90 C, 3 5 -dimer AG of -8.0 
kcal/mol, 3 '-terminal stability range of -4.8 to 11.6 kcal/mol, GC clamp stability of -8.0 kcal/mol, 
minimal acceptable loop AG of -1.9 kcal/mol, maximum number of acceptable sequence repeats of 6 and 
a maximum length of acceptable dimmers of 2 base pairs. In practice, when prepared, oligonucleotides 
are dried and re-hydrated in 3X SSC (a solution of sodium chloride and sodium citrate) at a concentration 
of 150 ng/|il. These particular values are the defaults for OLIGO 6 calculations. They can be adjusted by 
the user based on the physical biochemistry of the particular acids. These are good general values for 50- 
mers at this temp and CG content. 

[0063] In some embodiments, favorably evaluated target-length oligonucleotides 145, e.g., those found 
amenable to hybridization, are used as a query to a similarity search program 150, e.g., BLAST. The 
query is directed to a selected database 155 of genome data in order to determine whether the target- 
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length oligonucleotide is unique to the organism or organism set under investigation. To this end, 
preferred embodiments parse 150 the similarity search program output to identify oligonucletoides 
lacking significant similarity with other organisms in the selected database, e.g., unique target-length 
oligonucleotides 152. This is counter to the typical use of such search programs. For the purpose of this 
disclosure, "unique" or "uniqueness" is a function of thresholds, preferably controlled by the user, 
regarding identity, homology, score, expected (E) value and the length of the unique sequence under 
consideration. Identity, score, expected value are data returned in a typical BLAST search. In some 
embodiments, lacking significant similarity, e.g., "unique," means no BLAST hits or hits with a E-value 
less than le-5. In some embodiments, this determination of oligonucleotide uniqueness is conducted 
prior to evaluating the oligonucleotides for amenability to hybridization. 

[0064] At this point, the certain oligonucleotides that can be identified as unique to the set for which they 
were searched. Such oligonucleotides have utility in confirming the presence of at least one member of 
the set in both Bioinformatic and "wet" settings. Table 4 lists exemplary set-unique oligonucleotides for 
E. coli K12 identified by a method of this invention. Table 5 is an exemplary BLAST search result 
showing the sufficient uniqueness of an oligonucleotide of E. coli 0157:h7. These oligonucletides, with 
optional annotation, can be saved to a database C38 of unique genomic sequence, or otherwise added to 
the growing knowledge base of the genome of the organism under investigation. 

[0065] Referring again to Figure 1, unique oligonucleotide sequences were realigned 160 with the 
genome from which they were generated in order to determine the start and stop point of each sequence 
within the genome. In preferred embodiments, annotations within the genome in the region containing 
the unique sequence was transferred to the unique sequence. Annotated regions include sequences known 
to have a specific biological function such as protein coding regions, biologically active RNA encoding 
regions, promoter and regulatory elements, spacing elements within operons, protein binding sites, etc. 
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Oligo 
Number 


Organism . 


Sequence (5* to 3') 




E. coli K12 






NC 000913 29 part354 


91 




ACAGGATATAGTTATACCAGCGTTATTGTCGTTAGTGGTGAGTCGTCGGT 


92 




CAGTGATAATAACGTGACGCTGGATGGAAAGTTAACTGTTGTATCAGACA 


93 




TATAATCACAACGGTGAGCTTGTTA 1 1 ALvjCjCj 1 (jA 1 AA 1 CjL.(_ACA(jr 1 AAA 


94 




AGAGGTAGAAACCACCACAGGTGAAAAGAAAAATACCAATGTGACTTGTC 


95 




TATCACAGAGGTAGAAACCACCACAGGTGAAAAGAAAAATACCAATGTGA 


96 




TTATACCAGCGTTATTGTCGTTAGTGGTGAGTCGTCGGTATATCTGAATG 


Q 1 




TGAATATCACTGGTAACGTTCTGGTTGATAAGGATAAAACCGCAGACAAT 


98 




CAATACCGTTAATATGAATGGTGGACTTGAACTGATTGGAGAGAAAAACG 


99 




TTACGGCAGTGATAATAACGTGACGCTGGATGGAAAGTTAACTGTTGTAT 


100 




ATATTCGCCAGAAGGATTTATTATGCAAAGGAAAACTCTATTGTCGGCCT 


101 




G C AGATATC AC AGAGGTAGAAAC C ACCACAGGTGAAAAGAAAAATAC CAA 


102 




AGCGTTATTGTCGTTAGTGGTGAGTCGTCGGTATATCTGAATGGAGATAC 


103 




AAGAGCTTAAACGCTTACCCTCTGAATGCTCTCCTTTAGTCGAACAAAAC 


104 




GAAGTGGCAATATAATCACAACGGTGAGCTTGTTATTACGGGTGATAATG 


105 




CCGGCTTTCTTTATATTCGCCAGAAGGATTTATTATGCAAAGGAAAACTC 



TABLE 4 



BLASTN 2.2.7 [Jan-02-2004] 






RID: 1074620345-32204-105313520645.BLASTQ4 






Query= (50 letters) 






Database: All GenBank+EMBL+DDBJ+PDB sequences (but no EST, STS, GSS, or phase 0, 1 or 2 


HTGS sequences) 2,017,250 sequences; 9,771,1 19,756 total letters 






Distribution of 1 10 Blast Hits on the Query Sequence 








Score 


E 


Sequences producing significant alignments: 


(bits) 


Value 


gi|12519298|gb|AE005660.1|AE005660 Escherichia coli 0157:H7... 


100 


3e-19 


gi 1 1 3 3 64704|dbj | AP0025 69 . 1 1 Escherichia coli 0157:H7 DNA, c... 


100 


3e-19 


gi|262091|gb|S50878.1|S50878 tyrosine hydroxylase {5'region... 


40 


0.24 


gi|24430266|emb|AL928973.3| Mouse DNA sequence from clone R... 


38 


0.95 


gi|38083308|ref]XM_l 1 1684.3| Mus musculus similar to corneo... 


36 


3.8 


gi|38083112|reflXM 359287.1| Mus musculus similar to corneo... 


36 


3.8 


gi|33286907|gb|BCO55373.1| Mus musculus cDNA clone MGC:6076... 


36 


3.8 


gi|18412479|ref|NM_106620.1| Arabidopsis thaliana hypotheti... 


36 


3.8 


TABLE 5 
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[0066] Just as set-unique genomic sequences and oligonucleotides can be identified for strains of a 
particular organisms (E. coli as above), the method can be used to identify genomic sequences and 
oligonucleotides unique to higher level sets such as species, genus, family, clade, ad hoc sets, etc. 

[0067] Embodiments of the invention include both a method for profiling the hybridization response of 
set-unique target oligonucleotide sequences to a variety of organisms and sets of organisms and the 
resulting database of profiles. In these embodiments, target oligonucleotides are created in accordance 
with the unique sequences and are spotted to a microarray. Genomic DNA from each organism of interest 
is purified using procedures known to those skilled in the art. Preferably, the genomic DNA includes 
DNA from the area surrounding the region from which at least one target-length oligonucleotide was 
derived. In preferred embodiments, aliquots of the DNA are labeled, preferably with fluorescent dNTPs 
in a Klenow reaction. Each labeled DNA sample is allowed to hybridize with a microarray of the target- 
length oligonucleotides. The hybridized microarrays are wash, scanned, and the data is imported into a 
data visualization program (e.g., the suite of analysis software offered by Spotfire of Somerville MA). 
The data is evaluated to determine that target-length oligonucleotides exhibit true_positive and 
true_negative reactions to each organism of interest. 

[0068] In practice, multiple copies of target-length oligonucleotides and multiple trials of hybridization 
are used for each organism of interest. In this fashion, data can be collected on the variation of 
hybridization intensity to be expected between each target-length oligonucleotide and each organism of 
interest. For example, statistics regarding the distribution of hybridization intensity between a target- 
length oligonucleotide and its corresponding organism can be collected. Table 6 presents a notional 
example of hybridization intensity for two organisms against oligonucleotides found to be unique to 
separate strains of E. coli and the E. coli Shiga gene. Table 6 identifies exemplary unique oligonucleotide 
sequences and provides the hybridization intensity for those sequences showing the greatest differential in 
hybridization intensity as between E. coli K12 and E. coli ol57:h7 based on the query-length sequences 
NC_00913_29_part354 and NCJ)02695_194_part29 respectively, along with Shiga gene oligonucleotide 
that exhibit differential hybridization intensity between the two E. coli strains. Figure 3 illustrates this 
differential hybridization intensity between strains for oligonucleotides indicated by parenthetical 
numbers in Table 6. 
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Oligo 
Number 


Organism 


Sequence (5* to 3') 


Intensity 


r\A C7-LJ7 

Intensity 




E. coli K12, NC 000913 29 _part354 






(1) 91 




ACAGGATATAGTTATACCAGCGTTATTGTCGTTAGTGGTGAGTCGTCGGT 


33298 


13380 


92 




C AGTGATAATAACG TGACG CTGGATGGAAAGTT AACTGT TGT AT CAGACA 






(2) 93 




TATAATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAA 


139 to 


coon 

o8o0 


94 




AGAGGTAGAAACCACCACAGGTGAAAAGAAAAATACCAATGTGACTTGTC 






95 




TAT CAC AGAGGT AGAAACC AC C AC AGGTGAAAAGAAAAAT AC C AATGTGA 






(3) 96 




TTATACCAGCGTTATTGTCGTTAGTGGTGAGTCGTCGGTATATCTGAATG 


57432 


25699 


97 




TGAATATCACTGGTAACGTTCTGGTTGATAAGGATAAAACCGCAGACAAT 






(4) 98 




CAATACCGTTAATATGAATGGTGGACTTGAACTGATTGGAGAGAAAAACG 




bol4 


j 99 




TTACGGCAGTGATAATAACGTGACGCTGGATGGAAAGTTAACTGTTGTAT 






(5) 100 




ATATTCGCCAGAAGGATTTATTATGCAAAGGAAAACTCTATTGTCGGCCT 


1 t oo4 


04bz 


102 




AGCGTTATTGTCGTTAGTGGTGAGTCGTCGGTATATCTGAATGGAGATAC 






(6) 104 




GAAGTGGCAATATAATCACAACGGTGAGCTTGTTATTACGGGTGATAATG 


A\ OOOO 

1oobz 


2b4o 


[ 105 




CCGGCTTTCTTTATATTCGCCAGAAGGATTTATTATGCAAAGGAAAACTC 








E.coli 0157:h7, NC 002695 194_part29 






[ 278 










! 279 




ATTGTAAGGCGATTACTTTCTCAATCTTCTGAAATACGCAGTGCAAGTAG 






| 281 




TGTAAGGCGATTACTTTCTCAATCTTCTGAAATACGCAGTGCAAGTAGTC 






283 




TTAGGTCACACAATTCCAAAAGGCGGTTATGAACTTACTGAAGAACTTCC 






285 




CAATCTTCTGAAATACGCAGTGCAAGTAGTCTTCAGGTTGAAAGAATCTG 






(7) 287 




GTTATGAACTTACTGAAGAACTTCCCTGGGTTTTTGAAGCGATGAATGAG 


6038 


14850 


288 




GTCCAATTCTGTTAGGTCACACAATTCCAAAAGGCGGTTATGAACTTACT 






(8) 289 




GTAGTCTTCAGGTTGAAAGAATCTGGGTCGATGTTGATGACAGTTGGTAT 


14399 


47725 


292 




AAATCTCATTGTCGCTTCATGAAGGTCCAATTCTGTTAGGTCACACAATT 








E. coli Shic 


a Gene AF461172 






296 










297 




CTGACTATCATGGACAAGACTCTGTTCGTGTAGGAAGAATTTCTTTTGGA 






(9) 298 




AGGTACAACAGCGGTTACATTGTCTGGTGACAGTAGCTATACCACGTTAC 


9529 


18509 


(10) 299 




GTGAGCTATACGGAAAGTACACAAAAAGGAAGGTGCGACCACAATTAATA 


17781 


27191 


300 




ATCGCCATTCGTTGACTACTTCTTATCTGGATTTAATGTCGCATAGTGGA 






(11) 313 




CGGAAAGTACACAAAAAGGAAGGTGCGACCACAATTAATAACAAAATCTT 


7644 


20573 


314 




CGCTCTGCAATAGGTACTCCATTACAGACTATTTCATCAGGAGGTACGTC 






322 




ATTTACCAACAGATGGAATCTTCAGTCTCTTCTTCTCAGTGCGCAAATTA 






323 




AAGGAAGGTGCGACCACAATTAATAACAAAATCTTAAAAATTGCACATGG 






(12) 324 




AATTATTTACCAACAGATGGAATCTTCAGTCTCTTCTTCTCAGTGCGCAA 


5319 


20729 


(13) 325 




AAGTTATTTTTCGTTGACTCAGAATAGCTCAGTGAAAATAGCAGGCGGAG 


13167 


37404 



TABLE 6 



[0069] In one embodiment, the genomes of both E. coli K12 and E. coli 0157:H7 were investigated to 
identify unique sequence. This search included queries such as the query-length sequence in Table 1. 
Sequences such as NC_000913_29_part354 (Table l)and NC_002695_194_part29 and the E. coli Shiga 
Gene AF461 172 were BLAST searched against the entire NCBI nt database to verify that these sequences 
are unique as defined with respect to embodiments of the invention, e.g,, see Table 2 for the E. coli K12 
results. These results confirmed that the three query-length sequences noted directly above are unique. 
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The query length sequences were used to generate oligonucleotides. A exemplary list of these 
oligonucleotides is presented in Table 4 for E. coli K12. The full set of oligonucleotides were BLAST 
searched against the NCBI nr database to verify that they were unique at the oligonucleotide level; Table 
5 presents an exemplary extract of the BLAST output for an oligonucleotide of E. coli o 157:57. 
Oligonucleotides confirmed unique for the set were manufactured and spotted to microarrays. The 
microarrays were hybridized with genomic probes manufactured by Klenow labeling genomic DNA with 
cy3-dCTP. Oligonucleotides that demonstrated differential hybridization patterns were detected; see 
Table 6 for exemplary values and Figure 3 for a graphic representation. The thirteen oligonucleotides can 
now be used to distinguish between E. coli K12 and E. coli 0157:57. 

[0070] In some embodiments, the invention leverages the principles illustrated by the previous example. 
In those embodiments, the database of hybridization profiles is used as one source to pick 
oligonucleotides that are informative at decision points in decision tree. Structures appropriate for the 
decision tree include most taxonomic hierarchies, but any hierarchy where some oligonucleotides at 
sibling decision points have discernable differential hybridization will work. Placement of redundant 
unique and conserved phylogenetic specific oligonucleotide targets permits the identification of sets of 
organisms (e.g., family, class, order, genus, species, strain) specific branch points that can be used to 
identify organisms. In addition, placement of targets can include those sensitive to other features such as 
virulence genes, structural genes, biochemical genes, antibiotic resistance genes, housekeeping genes, 

[0071] Figure 4 illustrates a phylogenetic decision tree 400 in accordance with embodiments of the 
invention. In this particular embodiment, multiple oligonucleotide targets (including those from a 
database of unique oligonucleotides as described above that were found to provide useful differential 
hybridization between the organism sets under consideration) were spotted onto a microarray for each 
decision pint, e.g., 401 and exposed to a labeled complex sample for hybridization. Each decision point, 
e.g., 401 is scored with absent (-) or present (+) where: 




T 



A = number of (+) hybridized data points; 



B - number of (-) hybridized data points; and 
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[0072] In the embodiment illustrated in Figure 4, decision points having a score above a threshold are 
indicated as present (+). In other embodiments, actual hybridization values are reported. In further 
embodiments, absent/present calls are either made on the basis of statistics for all hybridization points 
associated with the decision point. In yet other embodiments, the node with highest score is called 
"present" under a branch point. 

[0073] In some embodiments, a confidence score is determined. In general each successive present call 
through the phylogenetic tree increases the confidence of the call below it and the confidence of the final 
identification. In Figure 4, the presence of two present calls for two different strains of E. coli decreases 
the confidence of an accurate identification at the strain level but increases the confidence of the 
identification at the level of the species. In some embodiments, the confidence score for a decision point 
is determined based on Bayesian analysis methodologies where present calls in the correct lineage 
contribute to the confidence score for the final identification. Factors that increase the confidence of the 
final identification include present calls in the phylogenetic positions in the tier above the absent/present 
call and additional absent calls laterally within that tier. Also contributing to the confidence of the final 
identification is the presence of lateral absent calls within the same tier as the absent/present call. Thus 
the absent/present call at any one location in the phylogenetic tree is dependent on all of the 
absent/present calls in all tiers. 

[0074] In some embodiments, the confidence score is represented as a numerical percentile between 0% 
and 100% by the following equation 

Confidence _ score = — xl 00 , where 

a = a + c + e , and 

p = b + d + f , where 

a = number of correct (+) calls in the correct lineage 
b = number of incorrect (-) calls in the correct lineage; 
c = number of correct (-) calls in the correct lineage; 
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d = number of incorrect (-) calls in the correct lineage; 
e = number of correct (-) calls in the incorrect lineage; and 
/ = number of incorrect (+) calls in the incorrect lineage. 

[0075] In some embodiments, a virulence score, a value used in conjunction with the final identification 
of the organism, is associated with the relative power and degree of pathogenicity of an organisms to 
produce disease. This value is derived from the number of virulence factors present, as detected by 
hybridization, compared to the total number of virulence factors represented in the array for a particular 
organism. Virulence targets can be spotted to the microarray independent of targets for absence/presence. 
This score is represented as a numerical value between -1 and +1 by the following equation: 



V A - the number of (+) hybridized virulence-specific data points, and 

V B = the number of (-) hybridized virulence-specific data points. 

V T - the total number of virulence-specific targets associated with the decision point. 

[0076] In addition, the ability to identify virulence versus structural components by oligonucleotide- 
specific hybridization permits the identification of recombinant organisms that contain structural 
components of one organism and the virulence components of a different organism. In some situations, a 
large number of data points may be required to identify all species and strains of micro organisms that 
might be found in a complex biological sample. 

[0077] Some methods of the present invention are conducted to mitigate the effect of "noise" by pre- 
screening. In accordance with these methods, a background sample of interest is obtained, and nucleic 
acid sequences in the sample are amplified using random amplification and combined with a microarray 
for hybridization. Signals from amplification products that hybridize with the microarray are recorded to 
be discounted in subsequent analysis. 

[0078] These methods are, for example, suited to customs applications. Customs officials at ports of 
entry including airports, harbors, and country borders can utilize prescreening to screen food samples for 
commonly occurring pathogens such as E. coli, Salmonella typhi, Hepatitis A virus and the like. In 
pathogen-free samples the level of hybridization observed to know pathogens on the array is minimal, the 



Virulence 
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hybridization profile of such pathogen-free samples is then used as a baseline level to subsequently 
identify contaminated samples. 
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